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Abstract 

An evolution partial differential equation for the surface of a non-wetting single-crystal film in an 
attractive substrate potential is derived and used to study the dynamics of a pinhole for the varying 
initial depth of a pinhole and the strengths of the potential and the surface energy anisotropy. The 
results of the simulations demonstrate how the corresponding parameters may lead to complete or 
partial dewetting of the film. Anisotropy of the surface energy, through faceting of the pinhole walls, 
is found to most drastically affect the time to film rupture. In particular, the similations support 
the conjecture that the strong anisotropy is capable of the complete suppression of dewetting even 
when the attractive substrate potential is strong. 

PACS numbers: 68.55.-a 
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I. INTRODUCTION 



Under typical operation conditions the thin solid film devices heat up to temperatures as 
high as 1000°, without melting. A high temperature activates the mass transport by surface 
diffusion. Thus the film morphology changes, which often results in harmful impacts on a 
device structural integrity. On the other hand, it is well known that thin liquid films (where 
the mass transport is by bulk hydrodynamic flows) are subject to dewetting caused by the 
long-range film-substrate interactions (which are also called wetting interactions). Dewetting 
results in film rupture and the agglomeration of liquid in droplets. The wetting interactions 
contribute noticably to the liquid film dynamics when the thickness of a film is of the order 
of few microns. Since this characteristic length scale is of the order of ten nanometers for 



solid films 



,l2|, and due to a high temperature activation barrier for surface diffusion, the 
dewetting of solid thin films is not as common. The structural integrity of a film must be 
maintained for a device to function properly. Thus it is important to understand the process 
of dewetting and find means to control it. 

One material system which exhibits solid film dewetting is the silicon-on-insulator. It 
comprises the ultrathin silicon film deposited on the amorphous Si02 substrate. The main 
feature of such system is the absence of stress at the film-substrate interface. The experi- 
ments with the sub-lOnm Si films at 800-900 °C demonstrate that dewetting starts at 
pinholes in the Si planar surface. The pinholes may exist prior to the annealing, or they form 
shortly after the temperature is raised. The critical depth of the pinhole which is necessary 
to initiate dewetting is presently unknown (apparently, not all pinholes dewet), nor is known 
the kinetics of the pinhole deepening and the pinhole shape. 

Due to anisotropy of the film surface energy it is natural to expect that the dynamics 
of a pinhole and the rupture process are qualitatively different from dewetting and rupture 
of the liquid surface having an isotropic surface tension. In particular, since the crystal 
surface tends to facet when its orientation falls in the unstable range, it is not a priori clear 
that dewetting, when it is initiated by substrate attraction of the initially slightly deformed 
free surface, not exhibiting unstable orientations, will proceed to rupture. Somewhat sim- 
ilar situation occurs in strongly anisotropic grain boundary grooving. Interestingly, it was 



demonstrated 



3| that such groove can be smooth and have the same shape as the isotropic 



groove, but the growth rate is reduced by a factor that depends on the degree of anisotropy. 



2 



This paper attempts to clarify the isssue by studying the dynamics of the pinhole with 
the varying initial depth, and for different strengths of the potential and the surface energy 
anisotropy. 

There have been several theoretical models published on the dynamics of thin solid films, 
where wetting interactions are accounted for these continuum models are primarily 

designed for the description of a self-assembly process in the course of the epitaxial crystal 
growth. In contrast to this paper, these models assume that the film wets the substrate, 
meaning that dfi w /dh\ h=const > 0, where \i w is the wetting chemical potential, and h is the 
film height. Also, dewetting of a solid monolayer was recently studied 8] using a solid-on- 
solid model and the Kinetic Monte-Carlo simulations. Dewetting of a monolayer proceeds 
through the nucleation and growth of holes in the film, which are caused by the discreteness 
of the underlying crystalline lattice and the substrate roughness. This is not the case for 
thicker films considered here. 

For solid films there is currently several choices of the wetting model, including the two- 
layer wetting model 0, If], El and the van der Waals potential [3, Q- In this study the 



former model is used due to its generality. It has been demonstrated in Refs. [H 



12| that 



the two-layer wetting model is consistent with ab initio calculations, but presently there is 
no experimental studies in which the wetting interaction potential has been measured. 

The continuum model of this paper is based on the well-established physico-mathematical 
framework which originates from the works of Mullins 3, Q] and Herring [jj], [lf| half a 
century ago. Mullins derived the first partial differential equations for the description of the 
surface morphology changes by surface diffusion (capillarity), evaporation- recondensation 
and bulk diffusion. Herring was first to systematically study the dependence of the equilib- 
rium crystal shape [l^] on the surface energy anisotropy. The experiments [3] demonstrate 
that the equilibrium shape can be partially or completely faceted, with the atomically sharp 
edges and corners and straight facets emerging when the anisotropy is strong, or with the 
smooth junctions and curved facets emerging when the anisotropy is weak. Other authors 
studied the dynamics of the crystal surface evolution with the strongly anisotropic surface 
energy and ways of regularizing the mathematical ill-posedness of the evolution problem (and 
the removal of the associated fine-scale instability) [sj], [3] - [^J. The preferred method of 
regularization is the addition of the small curvature-dependent term to the surface energy 
This term penalizes the sharp corners and makes them rounded on a 
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small length scale 



16] 



II. PROBLEM STATEMENT 

A two-dimensional film with the free one-dimensional surface z = h(x, t) is assumed, 
where h is the height of the film surface above the substrate. The surface evolves by surface 
diffusion. 

Under these assumptions, the governing equation for h(x, t) has the form 

where Q is the atomic volume, D the adatoms diffusivity, v the adatoms surface density, k 
the Boltzmann constant, T the absolute temperature, and \i = fi K + \x w the surface chemical 
potential. Here \i K is the regular contribution due to the surface mean curvature k [13l. 133]. 
Also (1+/?.^)~ 1//2 = cos 9, where 9 is the angle that the unit surface normal makes with the [01] 
crystalline direction, along which is the z-axis. (The x-axis is along the [10] direction.) Thus 
9 measures the orientation of the surface with respect to the underlying crystal structure. 
Note that here and below the subscripts x, t and 9 denote differentiation. 
The wetting chemical potential 

where, in the two-layer wetting model, 

y(h) = 1 U \9) + ( la - 7 (/) W) exp {-hjt) (3) 

is the thickness-dependent surface energy density of the film-substrate interface. Here y s = 
const, is the surface energy density of the substrate in the absence of the film, and i is the 
characteristic wetting length. Thus 7(h) — > y^ as h — >• 00, and ^(h) — > y s as h — > 0. 
y^\9) is the energy density of the film surface, assumed strongly anisotropic. The strong 
anisotropy possibly leads to faceting. (Below, when a reference to the surface energy density 
is made, the word "density" is omitted for brevity). 
y^(9) is taken in the form 

7 (/) (#)=7o(l + e 7 cos4#) + ^ 2 , (4) 
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where 70 is the mean value of the film surface energy in the absence of the substrate potential 
(equivalently, the surface energy of a very thick film), e 7 determines the degree of anisotropy, 
and 5 is the small non-negative regularization parameter having units of energy. The 5- 



term in Eq ; 
anisotropy 



J3|) makes the evolution equation ([T]) mathematically well-posed for strong 



29| - 32], 19|, [20|. (The anisotropy is weak when < e 7 < 1/15 and strong 



when e 7 > 1/15. 5 = in the former case. In the latter case the polar plot of ^'(9) has 
cusps at the orientations that are missing from the e quili brium Wulff shape and the surface 
stiffness 7^ + 7^ is negative at these orientations jisl Q. Thus the evolution equation 
is ill-posed unless regularized {29), Q.) The form PJ) assumes that the surface energy is 
maximum in the [01] direction. 

The curvature contribution to the chemical potential has the form 

fMlt = n( 7 w + > r $-H)K, (5) 

where 7(^(6') has the form (jlj) with the mean value 70 replaced by the effective value 7o 
which reflects the presence of the substrate: 
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7o = T-/ [7o + (is - 7o) exp (-z/£)\ dz. (6) 
fto Jo 

Here h = const, is the unperturbed film height. The if -term in Eq. (jSJ) is non-zero only 
for strong anisotropy, and it ensures the continuity of the chemical potential across corners 
on the faceted crystal shape. The form of this term is [3] 

Finally, k = -h xx {l + h 2 x )-^ 2 . 

To non-dimensionalize the problem, h is chosen as the length scale, and Hq/D as the 
time scale. Also, let r = £/h . Substitution of \i K and fi w in Eq. gives 

h = B^- (qPV _ APf + PW _ A P( 2 )) , (8) 

where h, x and t are non-dimensional, B = Q 2 i"y /(kTh l ), A = <V(7o^o) is the non- 
dimensional regularization parameter, and q, Pk and Pw ' are the following quantities. 

q = 1 + r - 1^ (1 -exp(-l/r)) > 0, (9) 
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P(2) = JL 

K " dx 



d 

— — (1 — 15e 7 cos 460 d cos 9; 

ox 



_ K 3 + cos() _ / cos() _j 



cos#; 



(10) 
(11) 



P 



3r _1 as^ (l + /i 2 ) + a 2 + 2a 4 /i 2 + aih^] 



(12) 



p(2) 



^=^[--'0.(i + *i) + 

^ {7h 2 xx h x - 2h xxx (l+h x ))]. 



(13) 



Note that Pi, 1} and Pi, Zj are proportional to exp (— h/r). In Eqs. (|T2|) and (|T3|) . ai = 
T — 1 — 17e 7 , a 4 = T — 1 + lle 7 and a 5 = T — 1 + 5e 7 /3. Thus Eq. ([8]) 



,(2) 



r - i 



«2 



contains five parameters, B, r, e 7 , A and T = 7^/70 (where A = when < e 7 < 1/15). 

From Eqs. (12)- (ED) it is easy to see that the requirement of zero wetting between the film 
and the substrate, dfi w /dh\h=h < 0, is equivalent to requiring negative a%, i.e. T < 1 + e 7 . 
Thus in what follows a± < 0. 

Eq. ([8]) can be written also in the widely used small-slope approximation (SSA), where 
\d/dx\ = e < 1: 

h t = B^- (gP« - AP« + P W _ AP( 2 )) , (14) 
pW = Ai/i^ + A 2 h 2 xx h x + A 3 h xxx h 2 x ; (15) 
P (2) = h ■ fl6l 



p(l) 



p(2) 



exp (—h/r) 



[a 2 h xx h x (l - 5fr 2 ) 



"^a, (ai + (2a 3 - 3ai) h 2 x + a 6 ^)] ; 



exp (—h/r) 



In Eq. (USD Ai = 15e 7 - 1, A 2 = 3 - 285e 7 and A 3 = 2 - 150e 7 . In Eq. (Jnj) a 3 



(17) 

(18) 
T - 1 + 3e^ 



and a§ = T — 1 — 25e 7 . Also note that the non- negative Ai signals that the surface energy 
anisotropy is strong. Eq. ( IT4T) has the same form as the evolution equation (26) of Ref. [5J. 

In the simulations reported below, the following values of the physical parameters are 
used: D = 1.5 x 1CT 6 cm 2 /s, fi = 2x 1CT 23 cm 3 , 70 = 10 3 erg/cm 2 , 7^ = 5 x 10 2 erg/cm 2 , 
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v = 10 15 cm -2 , kT = 1.12 x 1CT 13 erg, h = 10~ 6 cm, and 5 = 5 x 1CT 12 erg. These values 
translate into B = 3.57 x 10 -3 , T = 0.5 and A = 5 x 10~ 3 (or zero). Values of the parameters 
e 7 and r will be chosen later. 

The Method of Lines is adopted for the computation. Eqs. ([H]) and (TH|) are discretized by 



finite differences on a spatially-uni 



discretized in the nested fashion 



(1 2) 

'orm staggered grid in x. In Eq. ([S]) the terms Pk ' are 



34| . This, in particular, requires the explicit computation of 
only the first- and second-order partial derivatives h x and h xx . Thus, the numerical stiffness 
associated with the higher order derivatives becomes less pronounced, which manifests in the 
significant speed-up of the computation (in the cases of the large surface slope and faceting, 
the speed-up factor can be as large as 300). The results provided by this and by the standard 
scheme agree in all cases where such a comparison was made. The integration in time of 
the resulting coupled system of the ordinary differential equations is done using the implicit 
Runge-Kutta method. 

The initial condition in all runs is the surface 

h(x, 0) = 1 + a cos k max x, (19) 

where a is the non-dimensional perturbation amplitude (assumed not necessarily small, 
since the initial pinhole may extend deep into the film), and k max is the wavenumber for 
which the linear growth rate of the instability is the maximum. This wavenumber (or the 
corresponding wavelength \ max ) is usually referred to as "the most dangerous mode". In 
Eq. (fl9|) . < x < \ m ax (i-e., one wavelength). The boundary conditions at x = 0, \ ma x are 
periodic. 

Eq. (IS]) was used in all computations. Convergence was checked on finer grids. For small 
surface slopes the results provided by Eqs. ([HD and (TTJJ are identical. 

III. RESULTS 

The cases of weak and strong anisotropy are addressed separately. 

A. Zero or weak anisotropy 
Linearization of Eq. (fl4|) about the equilibrium solution h = 1 gives 

it = B (qA^ xxxx + exp (-l/r)r- 2 ai U) , Ai < 0, (20) 



where t) = exp (out + ikx) is the perturbation of the equilibrium solution, uj is the growth 
rate and k is the wavenumber. Thus for the growth rate one obtains 

oj(k) = B (qAik 4 - W k 2 ) , (21) 

where the notation is introduced Wq = exp (— l/r)r~ 2 ai. Obviously, Wq < for a\ < 0. It 
follows from Eq. (I2TI) that the equilibrium surface is unstable {uj{k) > 0) to perturbations 
with the wavenumbers < k < k c (where k c = {Wq/ (BqAi)) 1 ^ 2 ), co(0) = uj(k c ) = 
and, in the interval < k < k c , u(k) has the maximum value u max = — Wq /(ABqAi) at 
k ma x — k c /^/2. For k > k c the equilibrium surface is stable. Fig. 1 shows the sketch of uj{k). 

Simulations are performed for e 7 = and 1/18, and r = 0.1, 1. The graphs of the 
non-dimensionalized, reduced two-layer wetting chemical potential fi w are shown in Fig. 2. 
For r = 0.1, there is no substrate influence at the length scales larger than the unperturbed 
film thickness; for r = 1, the substrate- film interaction matters at all length scales less than 
ten times the unperturbed film thickness. For r = 1, the reduced potential is less strong, 
but more uniform across the film. Note that all figures in the paper show non-dimensional 
quantities. 

When the anisotropy is zero, the film dewets for both strengths of the wetting potential, 
as expected (other values of r, such as r = 0.01 and r = 10 also have been tried). Thus for 
the remainder of the paper the focus is on the realistic case of non-zero anisotropy, where a 
non-trivial evolution emerges. 

Fig. 3 shows the evolution of the surface shape for e 7 = 1/18 and r = 0.1. This was 
obtained using Eq. ([8]) (where A = 0), with the initial amplitude a = 0.01 in Eq. ffl9l . The 
shapes are smooth. Recalling the scaling, one can see that the 10 nm film completely dewets 
in approximately 45 seconds at 800 K. This is consistent with the order of magnitude time 
scale in the experiments 

HQ. 

Fig. 4(a) shows ln(l — h mm ) vs. time, where h m i n is the 
height of the surface at the tip of the pinhole. (Note that at t = 0, 1 — h min = a). Clearly, 
the linear theory fails to correctly predict the growth rate after t ~ 10 7 (6 s.). (Computing 
with the "small-slope" Eq. (fl4|) instead of Eq. ([8]) gives the same curve in Fig. 4(a), but 
the maximum slope of the pinhole wall at the rupture time is three times larger in the case 
ofEq. (tHJ).) 

Fig. 5 shows the evolution of the surface shape for r = 1. The transient shape approaches 
the equilibrium shape, which is different from the unstable trivial equilibrium h — 1, and 
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reaches this new equilibrium at t ~ 2 x 10 4 , i.e. in a small fraction of a second. Fig. 4(b) 
presents the growth rate data for this case. The equilibrium was confirmed by numerically 
solving the steady-state (boundary value) fourth-order problem for Eq. (I14p (using the 
shooting method). Most importantly, however, for all values of a in (0,0.99999] the initial 
deformation of the planar surface evolves to the shown equilibrium shape. Thus the model 
predicts the absence of dewetting in this case. 

B. Strong anisotropy 

Linearization of Eq. ([14"]) about the equilibrium solution h = 1 gives 

it = B {qA 1 ^ xxxx + A£ xxxxxx + WoU) , Ai > 0, W < 0. (22) 

The corresponding growth rate 

u(k) = B (qAik 4 - Ak 6 - W k 2 ) . (23) 

The equilibrium surface is unstable (u>(k) > 0) to perturbations with the wavenumbers 
< k < k c , and k 2 c = gA x /(2A) + (g 2 Af/(4A 2 ) - WyA) 1 / 2 . For k > k c the equilibrium 
surface is stable. The maximum growth rate and the corresponding wavenumber can be 
easily obtained. 

1. e 7 = 0.0697. 

For this value of e 7 (which is larger than the critical value 1/15 by only 4.5%) and 
r = 0.1 or r = 1 the evolution of the pinhole shape is qualitatively similar to the previous 
case e 7 = 1/18 and r = 1. For both values of r the equilibrium and the transient shapes are 
smooth. 

For r = 0.1 the equilibrium shape has the minimum at z = 0.92. This equilibrium state 
attracts all perturbations with the initial amplitudes a in (0,0.995). For a = 0.995 or 
larger, the pinhole tip approaches the substrate until it touches it, i.e. the pinhole dewets. 
By comparing this case to Fig. 3, one can see that the tiny supercritical anisotropy is almost 
successful in the suppression of dewetting. 

For r = 1, the equilibrium shape has the minimum at z = 0.95, and even for a = 0.99999 
there is no dewetting. 
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2. e 7 = 1/12. 



For r = 0.1, there emerges a smooth equilibrium shape with the minimum at z — 0.97, 
which attracts the initial perturbations with the amplitudes a in (0, 0.5]. The transient shape 
exhibits smooth curved "facets" and facet junctions. These facets gradually disappear as 
the equilibrium shape is approached. If a is in (0.5,0.99999], the transient shapes are truly 
faceted, and each shape evolves about the initial shape, becoming more faceted with time. 
Thus again, for e 7 = 1/12 and r = 0.1 there is no dewetting. Fig. 6 shows the initial and the 
final shapes for a = 0.8. Note that the tip of the pinhole is not faceted. In the dimensional 
time, the shape changes in Fig. 6 are practically instantaneous. After the low-energy 
orientations (0° and 90°) have been formed the time steps required to track the evolution 
become very small, and the computation is interrupted. (Of course, this is the drawback 
of the current model based on the evolution equation for the surface height; better model 
would employ the par ametric representation of the surface and will allow surface slopes of 
the 90° and larger [35j. I am presently working on such model.) 

For r = 1, the smooth equilibrium shape has the minimum at z = 0.965. This shape 
has a very narrow domain of attraction; already for a = 0.1, the clear tendency of the film 
to dewet can be seen in Fig. 7. The computation was interrupted after the last shape in 
Figs. 7(a) and 7(b) was dumped, since the time steps less than 10 -10 are needed to compute 
further. Such small time steps are warranted since the slope of the last shape near the tip 
is close to 90° and the evolution is very fast. These Figures also suggest that as a increases, 
the tip of the pinhole elongates while the rest of the pinhole shape changes little (is "left 
behind" ). In other words, the dewetting dynamics is confined to the tip of the pinhole. Thus, 
unless the facets form on the transient shape (Fig. 7(a)), a very narrow slit-like dewetting 
front is expected to extend to the substrate rapidly. If the facets do form, then the tendency 
to dewet will be partially or completely suppressed. 

IV. DISCUSSION 

The complex scenarios of a thin solid film dewetting under the action of the two-layer 
wetting potential have been computed using the nonlinear continuum model. These scenarios 
depend on the strength of the surface energy anisotropy, the characteristic wetting length 
and on the depth of the pre-existing pinhole defect in the film surface. 
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A film with the weak surface energy anisotropy dewets easily when the wetting potential 
is confined to the length scales of the order of the initial film thickness. When the wetting 
potential extends beyond the initial film thickness, there is no dewetting but instead the 
equilibrium emerges, where the pinhole is elongated but its tip does not touch the substrate. 

When the anisotropy is strong, it is generally expected that the faceting prevails over 
dewetting. For a very strong anisotropy either the smooth equilibrium emerges for shallow 
initial depths, or the transient shape becomes faceted and the dewetting ceases. Eqs. f fT2|) 
and ([13]) make clear the reason for such behavior. Indeed, the dominant terms in pffl (Pffl) 
are proportional to h~ l (h x 7 ). Thus for large slopes the "wetting" terms are small and they 
loose in the competition with the anisotropic terms and P^ ■ 

Two caveats must be noted. First, the averaging employed for the calculation of \i K 
is the approximation designed in order to obtain the simple model. Without averaging, 
the expression for the surface stiffness contains the exponential terms. As a result, the 
expressions for pjf* and, especially, Pk in Eq. (E]) are cumbersome and even the linear 
stability analysis for the strong anisotropy is complicated. In Appendix the partial proof 
is presented that use of averaging does not change the results qualitatively and that even 
the quantitative impact on the dynamics due to averaging is very small. Second, this model 
disregards the wetting stress However, as this reference suggests, the impact of such 
stress is not strong (except near the rupture point, where the stress singularity is expected). 
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Appendix 

Here Eqs. ([3]), (jlj) and ([51) are used as shown to derive the primary "wetting" part of \x K 
in the small slope approximation. The contribution of wetting energy to the regularization 
term, P\ in Eq. f fT6l) . is ignored. Thus in this extended SSA model Eqs. flT6|) - f|T8|) are 
unchanged, but P^ shown in Eq. (fl5l) is augmented by the exponential term 

P« = exp(-h/r)[-(A 1 +T)h xxx + 

(r" 1 (Ax + T) - A 2 h xx + 3Th xx ) h xx h x + 

(2T - A 3 ) h xxx h 2 x + r- 1 (A 3 - 2T) h xx hl] . (24) 
11 



Also, since averaging is not employed, q is set equal to one in Eq. f|T4|) . 

For the basic and the exdended models Fig. 8 compares the dewetting rates in the 
case e 7 = 1/18 and r = 0.1. This case was chosen because it is dominated by "wetting" 
contributions. The former rate is also shown in Fig. 4(a). Clearly the results are identical. 
For the large slope cases such as shown in Figs. 6 and 7, the differences are expected very 
small, if any, since the wetting correction to in Eq. ({TO]) is proportional to h~ 3 (vs. the 
h~ l dependence in Eq. ffT2l)). 
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FIGURES CAPTIONS 



FIG. 1. Sketch of the linear growth rate u(k). Perturbations with wavenumbers < k < k c 
are unstable and may grow nonlinearly until the film ruptures. 

FIG. 2. The reduced nondimensional wetting chemical potential \i w = — (l/r)(T — 
1) exp (—h/r). This formula is obtained when 7/(0) is taken isotropic and h x is taken zero 
in Eq. (j2J). T = 0.5. Solid line: r = 0.1; Dashed line: r = 1. 

FIG. 3. Surface shapes for r = 0.1 and weak anisotropy (e 7 = 1/18). 

FIG. 4. (a) Growth rate data for Fig. 3. Slope of the solid line equals the growth rate of the 
instability. For comparison, the dashed line is In a + uj max t, where 0J max is the linear growth 
rate, (b) Same as (a), but for Fig. 5. 

FIG. 5. Surface shapes for r = 1 and weak anisotropy (e 7 = 1/18). 

FIG. 6. Transient faceting of the surface for strong anisotropy. e 7 = 1/12, r = 0.1 and 
a = 0.8. 

FIG. 7. Transient surface shapes for e 7 = 1/12, r = 1. (a) a = 0.1. (b) a = 0.25. 

FIG. 8. Comparison of the dewetting rates for e 7 = 1/18 and r = 0.1. Solid line: basic 
SSA model with averaging, given by Eqs. (fT4j) - (|18p . Dashed line: extended SSA model as 
described in Appendix. 
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